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Abstract 

Recent experiments show that the isospin-violating mass sphtting of the 
B mesons is very small, but the best fits with a QCD sum rule analysis give 
a splitting of at least 1.0 MeV. The isospin-violating mass splittings of the 
charmed mesons, on the other hand, are in agreement with experiment. In 
this letter we show that the inclusion of 2"^'^ kind QED penguin diagrams can 
account for this discrepancy within the errors in the QCD sum rule method. 

PACS numbers: 12.38.Lg, 13.40.-f, 14.40.Jz 
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1. Introduction 

The development of a gauge-invariant method for QED corrections to composite 
systems^ has led to a consistent treatment within the QCD sum rule method of 
the three sources of isospin mass splittings: the current quark mass differences in 
the QCD Lagrangian, the nonperturbative QCD isospin violations which arise from 
u-d flavor dependence of vacuum condensates, and electromagnetic effects. The ap- 
plication of the method of Ref . [Q removed serious difficulties in earlier work on the 
isospin splittings of heavy-light mesonsp|. In a recent application^] to the charm 
and bottom pseodoscalar and vector mesons, satisfactory agreement with the exper- 
imental mass splittings was obtained within the expected accuracy of the method. 
However, although the D and D* isospin mass splittings were found consistent with 
experimental data, the theoretical result for the B''"-B° mass difference was about -1.2 
MeV, while the experimental value is 0.35 ± 0.29 MeVf^. Since the nonperturbative 
QCD effects tend to be quite small for the B mesons, this difference between theory 
and experiment seems to be somewhat larger than expected. 

Several years ago it was pointed out[P] that there is a novel electromagnetic ef- 
fect which modifies the quark-gluon vertex in a manner analogous to the penguin 
mechanism which leads to a weak correction to the quark-gluon vertex. This vertex 
modification is referred to in Ref. [^] as an "electromagnetic penguin of the second 
kind", because the term "electromagnetic penguin" was previously applied to the 
weak corrections to the quark-photon vertex. From the analytic form of the vertex 
modification 0, it is evident the the effect become increasingly important with in- 
creasing quark mass. We explore here the possibility that the neglect of this (second) 
QED penguin mechanism is the source of the discrepancy in the B isospin splitting 
found in Ref.[^. 
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The main results of the present paper are to point out that it is straightforward 
to find the largest nonperturbative QCD/QED effect as well as perturbative effects 
arising from this QED penguin mechanism. We estimate the contribution of nonper- 
turbative processes to the isospin- violating mass splittings of the D and B mesons and 
find a contribution of approximately 0.5 MeV to the pseudoscalar B-meson isospin 
mass splitting, while the penguins give a negligible contribution to the D mesons. 
The perturbative contributions are somewhat larger. The result is a 1.2-2.0 MeV 
mass splitting for the B and about 0.1 MeV for the D pseudoscalar mesons. From 
this we conclude that the QED penguins are of the correct magnitude to provide an 
explanation for the earlier discrepancy. 

2. Nonperturbative QCD/QED Penguin Isospin-Violating Mass 
Splitting 

Since the only isospin symmetry violating mechanism being considered in the present 
work is the (second) QED penguin vertex, it is convenient to take as the starting 
point of the QCD sum rules the correlators used for the heavy-light quark masses 
[without isospin violations] 0] . For the pseudoscalar current J5{x) =: q{x)i'j5Q{x) : 
the pseudoscalar current correlator is 

Mp) = ^ J dV^o\T{Mx)4{o))\o), (1) 

and 

n^.(p) = ^ / ciW^^(0|T(T/^(x)K(0))|0) (2) 

is the correlator for the vector current V^{x) =: g(x)7^Q(x) :. In the QCD sum rule 
method the correlator 11 (p^) is evaluated in two ways: 1) it is calculated starting from 



3 



QCD using the operator product expansion, and 2) it is treated phenomenologically 
by a dispersion relation. 

In the QCD calculation, the QED penguin isospin violations are obtained by 
including the quark-gluon vertex modification, depicted in Fig. 1, 

g"'^ = V,g, (3) 

at appropriate places in the QCD evaluation. In Eq. |^, g is the quark-gluon coupling 
constant and Vg, derived in Ref.0, is the penguin vertex modification. 

Keeping terms up to dimension D = 5, the microscopic nonperturbative QCD 
evaluation of n(p^) can be written as 

n(p2) = CjI + C,{qq) + C,{asG^) + C,{qi^-G)q) (4) 

where the coefficient C/ comes from the short distance correlation calculated using 
perturbative quark propagators, {qq) is the quark condensate, {agG"^) is the gluon 
condensate, (g(cr ■ G)q) is the mixed condensate and C„ are the Wilson coefficients. 
Cj is calculated from the process shown in Fig. 2a plus other two-loop diagrams. 
Since the graphs associated with the heavy-quark penguin vertices do not contribute 
to the mass splittings of interest in the present work, C/ can be obtained by a simple 
modification of previous calculations 0. Introducing the variable oj 



q^-M^ 



2Mq 

one can show that the coefficients C-i and Ca are of the form 



(5) 



Cs ~ -, (6) 

UJ 
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Since we eliminate these D=3 and D=4 processes in our sum rules, as we show below, 
we do not give the detailed results for C3 and C4. The coefficient C5 from the process 
of Fig. 2b is 



, . lim 
Bf[uj) = 



/(^), (9) 



1 

for the pseudoscalar currents, and vanishes for the vector current. For the process of 
Fig. 2c the D=5 mixed condensate coefficient is 

for both the pseudoscalar and vector currents. 
Defining the Borel transformation by 

/ d 
n! \ du 

and taking a derivative with respect to u;^ one obtains for the nonperturbative pen- 
guin graphs through D=5 

duB ~ 16^1 Mq ' ^'^^ 

for the pseudoscalar and 

duB 16uj% 

for the vector current. The Borel transformation of the D=3 and D=4 contributions 
[see Eq.|| are independent of ub and therefore do not contribute. Since we are only- 
considering isospin mass splittings in the present paper, the heavy quark vertex mod- 
ification, Vq, does not contribute, and we can use the zero-mass vertex modification 
from Ref.|p for Vg. The vertex modification for Vq in the limit of characteristic 
momenta larger than the Aqcd is0 

Vq = -^elil + ^ + lnb + lnb'), (12) 
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with h=-f^. (Note that Aqcd is simply the relevant hadronic scale independent of 

Q 

quark masses and is not necessarily directly related to standard QCD quantities such 
as Ajvj^.) The characteristic momentum of the heavy-light quark systems has been 
taken as the mass of the heavy quark, Mq. 

After the Borel transformation [Eq. the phenomenological forms of the corre- 
lators for the pseudoscalar and vector currents are 



and 

'phys ~ ^-"^M, 



n;hv,s = /.'77irre-^, (14) 



M'^—M'^ 

respectively, where A = , M is the mass of the meson states, and uo^ and uj^ 

are the pseudoscalar and vector Borel parameters. Defining II'^^-' as the correlator 
without the penguin processes, our sum rule for the mass shift due to the penguins, 
AM, is obtained by taking d(n - n(°))/da;B, 

for the pseudoscalar and 



V 



for the vector meson. The penguin expressions on the right-hand side of Eqs. |T5| and 



T6| consist of the sum of the D=5 nonperturbative terms given in Eqs. [T^ and |TT] and 



the perturpative contributions given by 



- ' dse ImCiis), (17) 



where uj is given by Eq. ^ with s=q^, and Sq is the threshold parameter for the 
continuum. The expression for the perturbative correlator, C/, can be obtained from 
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the results given in Ref0] and are found in Ref0. 

ry V Q 

ImCiis) = ^s{l-xff- + 2l{x)+ln{x)ln{{l-x) (18) 

5 1 5 

+ X - Y^)l'n{x) - {-- x)ln{l - x)], 

where x=Mq/s and l{x) is the Spence function||^. 

We use the results of the sum rule analysis without the penguins Q to fix the 
parameters needed for our estimate of the mass shifts due to the QED penguins. 
For the D meson: M^^ = 1.867, Mq = 1.5, Fd = 0.13 and ud = 0.7, all in GeV. 
For the B meson: = 5.279, Mq = 5.0, f^ = 0.095 and ujb = 0.63. Neglecting 
the perturbative contributions and keeping the nonperturbative contributions up to 
D=5, we obtain our results for the nonperturbative QED penquin mass shifts from 
Eqs.|lO|,|Tl,|l|,0 

[AM{D+) - AM(D°)]"P = OmMeV (19) 
[AM(S+) - AM(S°)]"P = 0.5MeV. 

The perturbative contributions are sensitive to the choice of the threshold parameter, 
Sq. The values of the threshold parameter which we use are taken from a study of 
the continuum in Refg]. The are 1.8-2.1 GeV^ for the D and 5.0-6.3 GeV^ for the B 
systems. We find for the perturbative contributions: 

[AM{D+) - AM(D°)]P* = 0.1 - 0.2Mey (20) 
[AM(S+) - AA'f (5°)]P* = 0.7 - l.bMeV 

for this range of the sq values. Note that the largest nonperturbative QED penguin 
contributions are negligible for the D isospin mass splittings, while they are of the 
order of the discrepancy in Ref. for the B isospin mass splittings. The net B mass 
splitting resulting from combining the nonpenguin result, M(B+)-M(B°)=-1.2 MeV 



from Ref. [H, with the penguin results of Eqs.|l^ and |2y is 

[AM(fi+) - AM(fi°)] = 0.0 - O.SMeV, (21) 
in agreement with experiment P] 
3. Conclusions 

In the present analysis we have shown that the QED penguins are of the proper 
magnitude to account for the experimental result that there is essentially no isospin 
mass splitting in the pseudoscalar B meson systems. From Eq.^, it is evident that 
the QED penquin contributions are much larger for the B than for the D systems, 
and that the good results of Ref . for the D system are not modified. We conclude 
that the QCD sum rule method can account for the heavy-light meson isospin mass 
splitting with the parameters that have been used for light quark systems. 

This work is supported in part by National Science Foundation grant PHY- 
9319641 and in part by the Department of Energy. 
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Figure Captions 

1. The QED penguin vertex. The oscillating curve represents a photon and the 

corkscrew curve represents a gluon. The straight line represents a quark. 

2. a) The perturbative QED penguin diagram; b), c) the D=5 nonperturbative 
penguin diagrams. The heavy line represents the massive quark (c or 6) and the light 
line represents the light quark {u or d). 
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